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Abstract 

We investigate the Generalized Parton Distributions (GPDs) of proton by expressing them in 
terms of overlaps of light front wave functions (LFWFs) using a simulated model which is able 
to qualitatively improve the convergence near the end points of x. We study the spin non-flip 
H{x,C,t) and spin flip E(x,C,t) part of GPDs for the particle conserving n ^ n overlap in the 
DGLAP region {(^ < x < 1). The Fourier transform (FT) of the GPDs w.r.t. to the transverse 
momentum transfer as well the FT of the GPDs w.r.t. C have also been obtained giving the 
distribution of partons in the transverse position space and the distribution in the longitudinal 
position space respectively. Diffraction pattern is obtained for both IL{x,cr,t) and £{x,a,t) in the 
longitudinal position space. 
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I. INTRODUCTION 


The Compton scattering of coherent light on an object is one of the most elementary 
processes of physics. In a general way, by measuring the angular and energy distributions 
of the scattered light information, the internal structure and the shape of the probed object 
can be accessed. With the advent of intense multi-GeV lepton beam facilities, it has become 
possible to experimentally study the Compton scattering at the smallest dimensions of mat¬ 
ter: the nucleon at quark and gluon level, where it is called Deep Virtual Compton scattering 
(DVCS). The term virtual here has the meaning that the incoming photon is radiated from 
a lepton beam, which presents the additional advantage of varying its 3-momentum inde¬ 
pendently of its energy. DVCS process 7 *(g) + p{P) —)■ +p{P'), where the virtuality 

of the initial photon = —Q^ is much large compared to the squared momentum transfer 
t = —{P — P'Y, provides a valuable probe to the structure of the proton. 

Within the framework of Quantum Chromodynamics (QCD), we intercept such reac¬ 
tion at the partonic level through the concept of Generalized Parton Distributions (GPDs) 
They are experimentally accessed through the overlap of DVCS and Beith-Heitler 


fl. 


Several experiments, such as, HI 


as exclusive vector meson production 

, I?!, ZEUS collaboration [sl, [oj and fixed target experiments at HERMES 


process as wel 
collaboration 6 

have hnished taking data on DVCS. Experiments are also being done at JLAB, Hall 
A and B 11] and COMPASS at CERN 12] to access GPDs. GPDs not only allow us to 
access partonic conhgurations with a given longitudinal momentum fraction, similar to deep 
inelastic scattering (DIS), but also at specific (transverse) location inside the hadron. GPDs 
depend on three variables x, Q t, where x is the fraction of momentum transferred, ^ gives 
the longitudinal momentum transfer and t is the square of the momentum transfer in the 
process. However, it has to be realized that only two of these variables ( (fully defined by 
detecting the scattered lepton ( = Xb, where Xb is the Bjorken variable used in DIS) and t 
(fully defined by detecting either the recoil proton or meson) are accessible experimentally. 
GPDs are much richer in content and carry much more information about the hadron struc¬ 
ture than ordinary parton distributions. In the forward limit of zero momentum transfer, 
the GPDs reduce to ordinary parton distributions. The GPDs give interesting information 
about the spin and orbital angular momentum of the quarks and gluons in the nucleon. 

It is well known that ( represents the longitudinal momentum transfer in the DVGS 
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process and in particnlar C = 0 represents the momentnm transfer only in the transverse 


direction. In Ref. 


13 


16| , GPDs have been stndied for C = 0 and it has also been shown that 


Fonrier transform of GPDs w.r.t Aj^ at (^ = 0 gives ipdpdfs (impact parameter dependent 
parton distribntion fnnctions). Since, in the DVGS process, the transformation of a virtnal 
photon into real photon reqnires a hnite transfer of longitndinal momentnm and the experi¬ 
ments always probe C 7 ^ 0, it becomes desirable to develop a deeper nnderstanding of GPDs 
at C 7 ^ 0. The sitnation for C 7 ^ 0 also has an important difference from C = 0 becanse as the 
proton loses longitndinal momentnm its transverse position is shifted by an amonnt which 


is proportional to C 


17l | . Farther, the transverse position of the proton is the vector snm of 


the transverse positions of its partons weighted by their longitndinal momentnm fractions. 
Therefore, when the proton loses its momentnm its transverse position get shifted from its 
original place. Even thongh the polynomiality conditions have to be satished by the GPDs 
with (C 7 ^ 0, it is very difficnlt to obtain a snitable parametrization of higher Fock states of 
the wave fnnctions in certain models. This is possible only if one considers the light front 
wave fnnctions (LFWFs) of simple spin 1/2 objects like dressed qnark or a dressed electron 
in theory. 

GPDs can be expressed as overlaps of LGWFs of the target hadron in light front gange 


18N20||. We have generalized the framework of QED by assigning a mass M to external 


electrons in the Gompton scattering process, bn 


; a different mass m to the internal electron 


line and a mass A to the internal photon line 2 l|]. The idea behind this is to model the 


strnctnre of a composite fermion state with a mass M by a fermion and a vector constitnent 
with respective masses m and A. These are off-forward overlaps in general and one reqnires 
not only particle conserving n ^ n overlap similar to forward pdfs bnt also n + 1 ^ n — 1 
overlap. Therefore, a spin 1/2 system can be represented as a composite of a spin 1/2 fermion 


and spin 1 vector boson with arbitrary masses 


22 


23l |. This one loop model is self consistent 


since it has correct correlation of different Fock components of the state as given by light front 
eigen valne eqnation. This model has been nsed to calculate the spin and orbital angular 
momentum of a composite relativistic system as well as the GPDs in impact parameter space 
and it also gives the Schwinger anomalous magnetic moment and the corresponding Dirac’s 
and Pauli form factor including the vanishing of the anomalous gravitomagnetic moment 
B(0). It provides a template for the wave function of an effective quark-diquark model of 
the valence Fock state of the proton LFWF. We would like to emphasize here that this is 
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extremely difficult to achieve in phenomenological models. In particular, two or three body 
Fock components can be obtained analytically from QED. Further, by taking the Fourier 
conjugate h_[_ (impact parameter) of the transverse momentum transfer Aj^, the GPDs can 
be expressed in impact parameter space. 

In Ref. [22, [2^ , the model with simulated bound state has been discussed. It is clear from 


these calculations that without involving the simulated bound state of the wave functions, 
one has to introduce the cut off on the momenta in the integrals for the expressions of GPDs. 
But differentiating the wave function w.r.t bound state mass improves the behavior of 
wave functions near the end points of x as well as at high k^. Differentiating the wave 
function w.r.t. produces a meson like behavior but it does not represent a model for 


meson. Even though some work has been done in the simulated model of LFWFs 


24l | where 


GPDs were discussed in impact parameter space i.e. by considering C = 0 but no work has 
been done so far for the case of C 7 ^ 0. There are three distinct regions in x. In the domain 
where ^ < x < 1 , there are diagonal 2 —)• 2 overlap contributions for both helicity non- 
flip H 2 ^ 2 {x,C,t) and helicity flip E 2 ^ 2 {xX^t). The GPDs H 2 ^ 2 {x,C^t) and E 2 ^ 2 {x,C^t) 
vanishes in the domain C — l<x<0. In the domain 0 < x < C, we have contribution from 
non diagonal n + 1 ^ n — 1 overlap. Since the differentiation of the single particle LFWF 
gives a vanishing result, this contribution vanishes in this model. It has also been shown that 
the DVGS amplitude expressed in terms of the variable a show diffraction pattern analogous 
to diffractive scattering of a wave in optics where the distribution in a measures the physical 


size of the scattering center in a 1 -D system 22, |23j, |25|] . The hnite size of the C, integration 


of the Fourier transform acts as a slit width and produces the diffraction pattern. Therefore, 
it becomes interesting to study the GPDs and the Fourier transform (FT) of the GPDs in 
the transverse as well as longitudinal position space. 

In the present work, we have calculated the El and E for non zero skewness. To understand 
the signihcance of the momentum fraction carried by the quarks in the process, we have 
considered the cases for hxed x and different values of C, as well as for hxed C, and different 
values of x. Both GPDs H and E are studied in longitudinal and transverse position 
space by taking Fourier transform w.r.t. C and Aj^ giving the distribution of partons in the 
longitudinal position space and the distribution in the transverse position space respectively. 
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II. KINEMATICS OF DVCS AND GENERALIZED PARTON DISTRIBUTIONS 


For the sake of completeness, we present the essential kinematics of deeply virtual Comp¬ 


ton scattering 


18l 


i*{q) + p{P) i{q') + p{P') • (1) 

The frame is specihed by choosing a convenient parametrization of the light-cone coordinates 
for the initial and hnal proton: 


P = P+, 0^, 


p+ 


P' = (1-C)P+, -Au, 


M2 + 




(1-C)P+^ 

where M is the proton mass. The four momentum transfer from the target is 

t + A2 


( 2 ) 

(3) 


A = P - F = CP+, A^, 




(4) 


where t = A^. In addition, overall energy-momentum conservation requires A =P —P' , 
which connects A^^, ^ and t as follows 

C^M^ + Ai 


t = 2P- A = 


1-C 


(5) 


The generalized form factors H, E are dehned through the matrix elements of the bilinear 
vector currents on the light cone |^: 


dy- 


Stt 


,^i.p*r/2 (F'l.v;(0)7+V-(9) |P) 




( 6 ) 


The off-forward matrix elements can be expressed as overlaps of the light front wave functions 


26| . For non-zero skewness ( there are diagonal parton number conserving contributions in 


the kinematical region (^ < a; < 1 and C — 1 < a; < 0 and there are parton number changing 
contributions in the region 0 < x < C- If we consider a spin 1/2 target state consisting of a 
spin 1 particle and a spin 1/2 particle, the contribution to the spin flip conserving and non 
conserving part of GPDs in the domain (/ < x < 1 can be expressed as 






2 - 1-2 






4(l-|)x/W 


p2^2(x,C,i) = 


5 














h' 'l./.t 




h' 'l./.t 




j Tj 


V'+i+i(a;', ^±)^|i+i(a;, ^±)+'0+i_i(a;', ^±)+^L|+i(a;', ^±)^Li+i(a;, ^±) 

(7) 




2M 


■E2^2{x,C,t) = 


(Pki_ r ,t* 


['*/'+1+1 ( 2 ;'> ^±)^|i+i(a^> + '0+i_i(a^'> 


t* 




where 


/ x-C r> r 7 -^a 

X = z -r, kj_ = k± — z -:A 1 . 


( 8 ) 


(9) 


1 -C’ " 1 -C 

We calculate the GPDs in simulated model of hadron LFWFs. To begin with, we present 
the two-particle wave function for spin up and spin down electron, which can be expressed 
as 


26 (: 


. 1 k_L) = -y/2{M - —)(p, 




X 


= 0 , 


( 10 ) 


and 


where 


V'+i+i(a;,fc±) = 0, 


-0+1 ^{x,k±) = -\/2{M -—)ip, 

+ 2 1 X 


'0ii+i(a^,^±) = -V^- 




-k^ + ik"^ 
{l-x) 

2 ^ x{l — X) 


(p{x,k±) = 


1 


vT 


X 


M2 - 


k^. +m2 fc2 _|__^2 


1—x 


( 11 ) 


( 12 ) 


As discussed earlier, we have used the generalized form of QED by assigning a mass M to 
the external electrons and a different mass m to the internal electron lines and a mass A to 
the internal photon lines. 
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FIG. 1: Two-particle LFWFs vs x for M = 150 MeV, m = 300 MeV, A = 300 MeV and fixed 
values of = fc in units of MeV. 


Differentiating the QED LFWFs with respect to improves the convergence of the 
wave functions at the end points of x as well as the behavior. We differentiate ip{x, k±) 
with respect to and have 


ip\x, k±) = 


VT 


X 


M2 - 


k^+m? k'^+X^ 

X 1—x 


(13) 


Using this, in FiglHwe have plotted the two particle LFWFs as a function of x for different 
values of k. It is clear from the plots that the wave functions converge at the end points of 
X and it’s behaviour is improved. 

The helicity non-flip GPD H can be calculated in this model, using Fqs. (7) and (8), 
and can be expressed as 


H{xX,t) 


1 — (/2 / x{x — C)(l — x) x^(l — x){x — 

( 1 - 0 * ( 1 - 0 * / 


A + -^2 + BI 3 + 
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E{x,C,t). (14) 


2{M{x — C) ~ ■'^(1 “ C)){Mx — m)x{l — x)^{x — () 

(1-0^ ' 


+ 




4(1-0 


The helicity flip term E can similarly be calculated from Eq. (|8]) and is given as 


E(x,C,t)=4M0l-C 


_(M - —)I^ ~ ~ 


77(1-0 {Mx-m) 

1 ( 

^M(x-C)-m(l-C)\\ 

VV 1 — X X 

1 — X ' 

V a:-C )) 


X 


(1-0^ 


(15) 


where 


^ f d‘^k± ^ (1 — Q;)dQ; 

^ i Jo Z 12 ’ 

f d‘^k± E ada 

= Iwr^’L 

r d‘^k± E a{l — a)da 

^ J IfZf ^""Jo Il3 ’ 

_ f d'^k±k± _ 1-x E y{l-yfA±dy 

^ “ J L\Li ~ I-c Jo g3 


(16) 


Li = k]_ — M x(l — x) + m (1 — x) + A x, 
L 2 


,2 M‘^ix — Q{l — x) „(1 —a^), , (1 — x)^A^ m^(l—x) 

k^ - ^ ^0. - - - 2- - -k^ ■ Ax + ^ + 


(1-0^ 


1-C 


(1-0^ (1-0 


A0(r-0 

1-C ’ 


D = 


(1 — x)'^A] 


a{l — a)— - ^ 0" + m^(l — a)(l — x) + A^(l — a)x — M^(l — a)x(l — x) 

(1-0^ 

M^a(x — (C)(l — x) ^ m^a(l — x) ^ A^Q!(x — C) 


( 1 - 0 ^ 


( 1-0 


1-C ’ 


Q = 


(1 - yfiX - xfAl ^ _ ^^^^ (1 -Q^^i _ M^{x -(){!- x) ^ m^(l -x) ^ 


B = 


1-C 

M^x(l — x) — 


(1-C)2 ^ -V (1-02 (1-C)2 

) + |/(m^(l — x) — M^x(l — x) + A^x), 


1-C 


^^(^“C)5 , , 2 


2 n 5 (1-O^l^i 5 x 2 , M2(x-C)(l-a;) m2(l-x) 


(1-0^ 


— m (1 — x) — A x + 


( 1 - 0 ^ 


1-C 


AQx-O 

1-C 


(17) 


To understand the behaviour of the GPDs with the change in momentum transfer, in Fig. 
[21 we have plotted the helicity non-flip GPD i7(x, (C, t) and helicity flip GPD i7(x, (C, t) as a 
function of x with fixed value of C, whereas in Fig. [3] we have plotted the helicity non-flip 
GPD i7(x, (C, t) and helicity flip GPD E{x, t) as a function of ( with hxed value of x. For 







































FIG. 2: Plot of i/ vs X for a fixed value of C with different values of —t and plot of FI vs x for a 
fixed value of C with different values of —t. Parameter t is in MeV"^ 




FIG. 3: Plot of H vs C for a fixed value of x with different values of —t and plot of FI vs C for a 
fixed value of x with different values of —t. Parameter t is in MeV‘^ 


the numerical calculations we have used m = 300 MeV, M = 150 MeV, and A = 300 MeV 
22l-l24| . A cursory look at the plots reveal that the behaviour of the GPDs is the same and 
is independent of the value of momentum transfer |f|. In Fig. [2l H{xX,t) and E{xX,t) 
increase with x, reach maximum and then decrease. It is important to mention here that 
since x is the momentum fraction of the active quark, at x = 1, the active quark carries 
all the momentum and the contribution from other partons is expected to be zero at this 
limit. It is also observed that the peak of H{xX,t) shifts towards higher value of x as |f| 
increases suggesting that the active quark is more likely to have a large momentum fraction. 


9 















































We are considering here only the leading Fock space component of the target hadron state. 
The higher Fock space components are more likely to contribnte in the small x region. The 
helicity flip GPD E{x, t), related with the qnark orbital angnlar momentnm, also decreases 
as |t| increases. In the case, the peak occnrs at a lower valne of x as compared to the peak in 
the case of H{x, (, t). However, the shifting of peak towards higher valne of x as |t| increases 
is similar to that in the case of H{x, (, t). 

On the other hand, in Fig. [3l when H{x,(,t) and E{x,(,t) are plotted as a fnnction 
of C with hxed valne of the fraction of momentnm carried by the active qnark x, the peak 
occnrs at same valne of (. However, the amplitnde decreases as the valne of |f| increases. 
Since ( represents the longitndinal momentnm transfer in the process and we are working in 
the region ^ < a; < 1, therefore, H{xX,t) increases as the valne of the momentnm transfer 
( increases and after reaching a maxima it starts decreasing. It is observed that even if 
we increase |t|, the peak still occnrs at same valne of (. This is dne to the fact that the 
actnal momentnm dnring the process is carried by the qnark and since it is fixed in this case, 
therefore any increase in ( simply implies the increase in longitndinal momentnm fraction. 
This does not affect the position of the peak which is observed at the same valne of ( even 
for different valnes of t. 

III. GPDS IN POSITION SPACE 

The FT with respect to the transverse momentnm transfer Aj_ gives the GPDs in trans¬ 
verse impact parameter space. We introdnce b± conjngate to Aj^ giving 



(18) 



(19) 


where A=|A_l|. Here b=\b±\ is the impact parameter measnring the transverse distance 
between the strnck parton and the center of momentnm of the hadron and is dehned as 
Y^Xibi = 0 (the snm being over the nnmber of partons). Since we are working in the limit 
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of non-zero skewness, the transverse location of the proton before and after scattering is 
different and this transverse shift depends upon ( and b±. Even when integrated over x the 
information does not wash away. In the DGLAP region ((^ < x < 1), the parameter b = |6x| 
gives the location of the quark where it is pulled out and put back to the nucleon whereas 
in ERBL region it describes the location of quark-antiquark pair inside the nucleon. In Fig. 
m we plot 'H(x,C,6_l) and £^(x,C,6_l) as a function of |6 _l|. We have hxed the value of x 
and have taken different values of C,. We know that C can never be zero in the experiments 
because their is always a hnite momentum transfer in the longitudinal direction. Therefore, 
probability interpretation is not possible. It is clear from the plots that the peak of GPDs 
in transverse impact parameter space decreases as value of C, increases. For the sake of 
completeness, in Fig. O we have plotted 'H{xX,b_L) and £{x,(,bi_) as a function of |6j_|. 
In this case however, we have hxed the value of ( and have taken different values of x. 
In contrast to Fig. IH in Fig. [5l the peak of GPDs in transverse impact parameter space 
increases x increases. The smearing in the b± space is due to the multiparticle correlation. 
As we increase the value of ( for a hxed value of x, smearing in the b± space increases 
implying that the momentum transfer in the longitudinal direction increases. 

Since we are ehectively looking at the initial and hnal waves of the scattered proton where 
the hnal state proton wave function gets modihed to the incident proton wave function 
because of the momentum transferred to the quark in Gompton scattering. The change in 
quark momentum along the longitudinal direction can be Fourier transformed to a shift in 
the light front position of the struck quark. Similar to the transverse momentum transfer. 


following Ref. 


22 


231], we introduce a longitudinal boost invariant impact parameter a, 


which is conjugate to longitudinal momentum transfer (. One can thus simulate a change 
in the quark’s longitudinal light front co-ordinate by an amount a = ^b~P^, designated as 
the scaling parameter. This is analogous to dihractive scattering of a wave in optics where 
a plays the role of the physical size of the scattering center in a one-dimensional system. 
We are using t to set the change in the transverse momentum of the quark in the scattering 
the scale used is dynamical in nature. The GPDs in longitudinal position space are now 
expressed as 

1 r 


^{x, a, t) = 


27r 


1 rG 




( 20 ) 

( 21 ) 
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FIG. 4: Fourier spectrum of 'H{x,Q,h±) vs |6 _l| for fixed value of x and different values of Q and 
Fourier spectrum of ^(x, vs |6 _l| for fixed value of x and different values of C- 




FIG. 5: Fourier spectrum of 'H{x,Q,h_\_) vs |6 _l| for fixed value of C and different values of x and 
Fonrier spectrum of vs |6 _l| for fixed value of C and different values of x. 


In the region ^ < x < 1, the npper limit of integration is given by (max 


(max 


2M2^\ 


1 + 


4M2 


1 ), 


( 22 ) 


whereas in the region 0 < x < C, the npper limit is given by the valne of x. 

In Fig. [6l we plot the Fonrier pattern for 'H(x, a, t) and S{x, a, t) as a fnnction of a in the 
longitudinal position space. Diffraction pattern is obtained for both ^{{x, a, t) and £{x, a, t) 
where (max plays the role of the slit width. Since the position of the minima (measured from 
the center of the diffraction pattern) are inversely proportional to the slit width, the minima 
move away from the center as the slit width ((max) decreases. In both the plots the primary 
maxima is followed by secondary maxima and as the value of |t| increases, the minima moves 
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FIG. 6: Fourier spectrum of T-L{x, a, t) vs a for fixed value x with different values of —t and Fourier 
spectrum of £{x,a,t) vs a for fixed value x with different values of —t. Parameter t is in MeV'^ 


towards smaller value of a. Therefore, the longitudinal size of parton distribution becomes 
longer and the shape of the conjugate light cone momentum distribution becomes narrower 
with increasing |t|. It can also be added that the position of minima is independent of 
helicity and it appears at same points in both cases. This is a characteristic of diffraction 
pattern obtained in single slit experiment. 

IV. CONCLUSIONS 

In the present work we have studied the GPDs in transverse and longitudinal position 
space. We have used a two body model which represents a composite system consisting of 
a fermion and vector boson with arbitrary masses. The calculations have been done with 
non-zero ( in the DGLAP region x > (. We have considered the n ^ n parton number 
conserving overlap. Differentiating the denominator of the light front wave function w.r.t 
improves the convergence near the end points of x and an improved behaviour is observed. 

The behaviour of the helicity non-flip GPD H{xX,t) and helicity flip GPD E{xX,t) as 
a function of x with hxed value of ( increases with x, reaches maximum and then decreases. 
The active quark carries all the momentum in the limit x —)■ 1 and the contribution from 
other partons is expected to be zero. The active quark is more likely to have a large 
momentum fraction as the value of |t| increases. We are considering here only the leading 
Fock space component of the target hadron state. The higher Fock space components are 
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more likely to contribute in the small x region. On the other hand, when H{xX,t) and 
E{xX, t) are plotted as a function of ( with hxed value of the fraction of momentum carried 
by the active quark x, the peak occurs at same value of (. This is due to the fact that the 
actual momentum during the process is carried by the quark and since it is hxed in this case, 
therefore any increase in ( simply implies the increase in longitudinal momentum fraction. 
This does not affect the position of the peak which is observed at the same value of ( even 
for different values of t. 

The impact parameter space representations are obtained by taking the Fourier transform 
of the GPDs with respect to transverse momentum transfer and it probes the parton distri¬ 
butions in impact parameter space. If |6j_| is the impact parameter measuring the transverse 
distance between the struck parton and the center of momentum of the hadron, the trans¬ 
verse location of the proton before and after scattering is different and this transverse shift 
depends upon ( and b±. The smearing in the b± space occurs due to the multiparticle cor¬ 
relation and increases with the increase of ( implying that the momentum transfer in the 
longitudinal direction increases. 

Further, if we introduce a longitudinal boost invariant impact parameter a conjugate to 
longitudinal momentum transfer after Fourier transform, we get the GPDs in longitudinal 
position space. Both H and E show the diffraction pattern analogous to diffractive scattering 
of a wave in optics where the distribution is in a space. Similar pattern has been observed in 
some other models also but the general features of this pattern are independent of specihc 
models and depends mainly on the hniteness of C, integration as well as the dependence 
of GPDs on X, C, and t. However, in order to get the full Lorentz invariant picture in 
longitudinal space, one has to study the GPDs in a; < C region as well. 
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